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The observed stability of nanobubbles contradicts the well-known result in classical nucleation the¬ 
ory, that the critical radius is both microscopic and thermodynamically unstable. Here nanoscopic 
stability is shown to be the combined result of two non-classical mechanisms. It is shown that the 
surface tension decreases with increasing supersaturation, and that this gives a nanoscopic critical 
radius. Whilst neither a free spherical bubble nor a hemispherical bubble mobile on an hydrophobic 
surface are stable, it is shown that an immobilized hemispherical bubble with a pinned contact rim 
is stable and that the total entropy is a maximum at the critical radius. 


I. INTRODUCTION: NANOBUBBLE TROUBLE 

Experimental evidence for nanobubbles was first pub¬ 
lished in 1994*i and since then their existence has 
been confirmed from various features of the mea¬ 
sured forces between hydrophobic surfaces,— ~- 4 includ¬ 
ing a reduced attraction in de-aerated water^r— and 
from images obtained with tapping mode atomic force 
microscopy— &iiL12r2L For a recent review of theory and 
experiment, see Ref. m 

In the fields of nucleation science and bubble research, 
three concepts are import ant d 9 : 20 First, the critical ra¬ 
dius, which is the point at which the derivatives of the to¬ 
tal entropy (equivalently the free energy) vanish. Second, 
the Laplace-Young equation, which says that the inter¬ 
nal gas pressure of the bubble is larger than the pressure 
of the surrounding liquid by an amount equal to twice 
the surface tension divided by radius. Third, for the case 
of heterogenous nucleation and bubbles adhered to sur¬ 
faces, the Young equation, which says that the cosine 
of the contact angle is equal to the difference in surface 
energies of the solid divided by the liquid-vapor surface 
tension. 

The trouble with nanobubbles is that 

1. the critical radius is an unstable equilibrium point 
(it is a saddle point), whereas nanobubbles appear 
to be stable 

2. the conventional critical radius for realistic param¬ 
eters is an order of magnitude or so larger than the 
typical nanobubble radius of curvature 

3. the Laplace-Young equation predicts that typical 
nanobubbles have an internal gas pressure one or 
two orders of magnitude larger than atmospheric 
pressure 

4. the measured contact angle for typical nanobubbles 
is significantly larger than the contact angle mea¬ 
sured for macroscopic droplets on the same surface 

Any explanation for the existence of nanobubbles must 
reconcile these four points with conventional thermody¬ 
namic principles and realistic physical models. 

The first point is such a strong qualitative contradic¬ 
tion with the experimentally observed stability (hours 
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FIG. 1: Total entropy of a spherical air bubble in water as a 
function of radius for a supersaturation ratio of 2 and a surface 
tension of 0.072 N/m. The solid curve follows the mechanical 
equilibrium path, dStot/dR = 0, the dashed curve follows 
the diffusive equilibrium path, dStot/dN = 0, and the dash- 
dotted curve follows the constant number path, N = N CI it. 
The dotted line marks the critical radius. 


or days) of nanobubbles that it calls into question con¬ 
ventional nucleation thermodynamics itself. The second 
point represents a quantitative discrepancy between the 
length scales for bubbles predicted by thermodynamics 
using conventional parameter values and the measured 
length scales of nanobubbles. The third point says that 
there is a strong thermodynamic gradient between the 
gas inside the nanobubble and the atmosphere, which 
should lead to the runaway dissolution of a nanobubble 
on time-scales orders of magnitude shorter than the ex¬ 
perimentally observed lifetimes. The fourth point again 
appears to prove that conventional macroscopic ther¬ 
modynamics, both bulk and surface, is inapplicable on 
nanoscopic length scales, which broad point is at vari¬ 
ance with diverse experimental measurements in other 
systems, molecular-level computer simulation data, and 
conventional wisdom. These four points arc unlikely to 
be independent of each other. 

Figure[l]shows the total entropy for a nucleating spher¬ 
ical bubble. A stable equilibrium point corresponds to a 
local maximum of the entropy. (In this article all thermo¬ 
dynamic results are cast in terms of entropy. These can 
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be translated into free energy by multiplying throughout 
by minus the temperature. Hence maximum total en¬ 
tropy is the same as minimum free energy.) It can be 
seen that when the growth of the bubble occurs along a 
path of either mechanical or diffusive equilibrium, the to¬ 
tal entropy is a minimum at the critical radius. However 
if the change in size occurs at either constant number 
(shown) or constant radius (not shown), then the total 
entropy is a maximum at the critical radius. Hence the 
total entropy has a saddle point at the critical radius, 
which means that it is a point of unstable equilibrium. 
The critical radius for these parameters is about an order 
of magnitude larger than the radius of curvature typically 
measured for nanobubbles. 

The physical reason for the instability of the critical 
radius is readily understood. At the critical radius the 
diffusive (number) and mechanical (radius) gradients are 
zero. A fluctuation to a larger radius along the path of 
mechanical equilibrium (mechanical equilibrium occurs 
faster than diffusive equilibrium with the reservoir be¬ 
yond the immediate vicinity of the surface) reduces the 
internal pressure and hence the gas chemical potential 
inside the bubble. There is now a thermodynamic gradi¬ 
ent that forces dissolved gas from the reservoir into the 
bubble, which increases the internal pressure. This cre¬ 
ates a mechanical force imbalance, which acts to increase 
the volume of the bubble at the expense of the volume 
of the reservoir. This increases the radius, magnifying 
the original fluctuation. The opposite occurs if the ini¬ 
tial fluctuation decreases the radius, in which case the 
bubble is driven to zero radius. 

An apparently trivial but essential point is easily over¬ 
looked in this explanation: The instability at the critical 
radius relies upon the fact that the volume of the bubble 
increases with increasing radius. It will turn out that this 
simple geometric fact holds the key to understanding the 
thermodynamic stability of nanobubbles. 

The purpose of this paper is to show what it takes 
to make the critical radius of a bubble both stable and 
nanoscopic. 


II. CONVENTIONAL NUCLEATION 
THERMODYNAMICS 

Consider an air bubble in water. Only the air need 
be explicitly accounted for, as the water solvent can be 
considered always at equilibrium (see Appendix El- The 
total entropy is taken to be the sum of three terms: the 
bulk entropy of the sub-system, which is the bubble, the 
bulk entropy of the reservoir, which is the air dissolved 
in water, and the interfacial entropy of their boundary^ 


S tot (N,V,T) 


7 A 


= S s (N,V,T) + S t (N t ,V t ,T) - 

pV pN Att^R 2 


= S a {N,V,T) - + 


T 


The total system is in thermal equilibrium at tempera¬ 
ture T. Here N is the number of molecules in the bubble, 
V = AnR 3 /3 is its volume, A = AnR 2 is its area, and R is 
its radius. The pressure of the reservoir is p, which equals 
atmospheric pressure, the chemical potential of the reser¬ 
voir is /i, and the liquid-vapor surface tension is 7 . As 
shown in Appendix E it would be more precise to set 
the external reservoir pressure equal to the partial vapor 
pressure of air at saturation, p = p' ss 0 . 8 p a t m , but for 
simplicity this is not done here. 

In the capillarity approximation, the bubble is taken 
to be an ideal ga o 19 i 21 
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where k B is Boltzmann’s constant and A is the thermal 
wave length. Although the entropy of a real gas differs 
from this quantitatively, the ideal gas model is not ex¬ 
pected to affect the results qualitatively. 

A bubble can only be in equilibrium if the water is su¬ 
persaturated with air (see below). The supersaturation 
ratio, s, is the ratio of the actual concentration of dis¬ 
solved air to the saturation concentration, which equals 
the ratio of the air pressures that would give the two 
concentrations, s = p as /p\ where the dagger here and 
below denotes saturation. The chemical potential can be 
obtained from the supersaturation ratio by the standard 
ideal gas result^ 


A 3 

A 4 = k B T In (2.3) 

k B l 

In this work the saturation air pressure will be taken 
to be atmospheric pressure, p' = p. This neglects the 
vapor pressure of water, which is about one fifth of an 
atmosphere. 

The number derivative is 


dStot{N, V, T) _ p s u 
dN T + T , 


(2.4) 


which vanishes when p s = p. Using the ideal gas equa¬ 
tion this is 


p = A -3 e »/k B T = JP_ 
H k B T’ 


(2.5) 


where p = N/V is the number density of air inside the 
bubble. This is the equation for diffusive equilibrium 
between the air inside the bubble and the air dissolved 
in the water. 

The volume derivative is 


dS tot {N,V,T) _ Ps p 2 7 

dV T T TR' 

This vanishes when 


( 2 . 6 ) 
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which is the Laplace-Young equation. This is the equa¬ 
tion for mechanical equilibrium for the bubble. 

When the bubble is in mechanical equilibrium, the 
Laplace-Young equation says that the internal pressure 
exceeds the external pressure. Since the pressure is a 
monotonic function of the chemical potential, and since 
the external pressure is equal to the saturation vapor 
pressure, the chemical potential of the air inside the bub¬ 
ble must exceed the saturation chemical potential. This 
is the reason that a bubble in mechanical equilibrium can 
be in diffusive equilibrium only with a supersaturated 
solution^ 

The critical radius is the point at which diffusive and 
mechanical equilibrium simultaneously hold. Using the 
the ideal gas equation of state, p s = pk^T, and substitut¬ 
ing the diffusive equilibrium equation into the mechanical 
equilibrium equation, the critical radius is 


R 


crit — 


27 

0- i)p' 


( 2 . 8 ) 


The critical radius scales linearly with the surface tension 
and decreases with increasing supersaturation ratio. The 
depth of the entropy minimum (equivalently, height of 
the free energy barrier) scales with the cube of the critical 

radius, *Stot,crit = l)P^f / crit/2. 

The numerical results in Fig. [T] show that the criti¬ 
cal radius is a point of unstable equilibrium, which is 
of course well-known in classical nucleation theory ] 19 ! 20 
See Appendix [B] for a mathematical proof of the re¬ 
sult. Also this, the conventional formula for the criti¬ 
cal radius, would require a supersaturation ratio s = 10 - 
100 (equivalently, water in equilibrium with air at 10 - 
100 atmospheres) in order to produce a critical radius 
-Rcrit = 10- 7 10 “ 8 m, which are the radii measured for 
nanobubbles. 


III. STATE DEPENDENT SURFACE TENSION 

As explained above, if a bubble is in diffusive equilib¬ 
rium with the water phase, then the solution must be 
supersaturated with aii > 22 The effects of supersaturation 
on the surface tension are neglected in the classical the¬ 
ory, but in fact these effects can be quite dramatic. 

The barrier to bubble growth is the depth of the en¬ 
tropy at the critical radius, which scales with the cube 
of the surface tension. At the spinodal supersaturation, 
there is no barrier to bubble growth, which means that 
the surface tension must vanish J2 Hence the surface ten¬ 
sion must decrease with increasing supersaturation* 22 r— 
This general thermodynamic result is consistent with 
what is known experimentally from high pressure mea¬ 
surements, namely that the surface tension of water is 
lower in the presence of oxygen or nitrogen and water 
vapor than in the presence of pure water vapor ^ 

Computer simulations of the supersaturated liquid- 
vapor interface show that to a reasonable approximation 
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FIG. 2: Total entropy of a spherical air bubble in water 
as a function of radius following the mechanical equilib¬ 
rium path. The solid curves are for a supersaturation ra¬ 
tio of s = 2 and the dashed curves are for s = 3. The 
lower curves of each pair use the saturation surface tension 
Y = 0.072 N/m, and the upper curves of each pair use the 
supersaturated surface tension (spinodal supersaturation ra¬ 
tio s* = 5), 7 ( 2 ) = 0.054 N/m, and 7 ( 3 ) = 0.036 N/m. The 
arrows mark the respective critical radii (local entropy mini¬ 
mum). 


the decrease is linear^ — 


7 ( s ) 



' t t' 
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Here 7 1 = 0.072 N/m is the saturation or coexistence sur¬ 
face tension (it was denoted plain 7 above), s = p s /p' g 
is the reservoir supersaturation ratio, which is the ratio 
of the supersaturation air vapor pressure to the satu¬ 
ration air pressure, and the dagger and double dagger 
denote the coexistence and the spinodal values, respec¬ 
tively. Obviously, s* > = 1 and > P g - Simulations 

of a Lennard-Jones fluid at various temperatures found 
that the surface tension was well-fitted by this with spin¬ 
odal supersaturation ratios in the range s* = 3-6 j22t— 

This surface tension may be called the supersaturated 
surface tension. It has a constant value independent of 
the radius of the bubble and dependent only on the su¬ 
persaturation ratio of the reservoir, (the concentration 
of air dissolved in the water divided by the saturation 
concentration). 

If one uses this reservoir state-dependent surface ten¬ 
sion, then all of the results of classical nucleation theory 
hold, with the value of the surface tension being replaced 
by the supersaturated value, 7 1 => y(s). This means that 
both the critical radius, which scales linearly with the 
surface tension, and the nucleation barrier, which scales 
cubicly with the surface tension, are reduced from their 
classical values. 

Figure [2] shows the effect of the supersaturated sur¬ 
face tension on the total entropy of a bubble. It can be 
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FIG. 3: Critical radius (log scale) as a function of supersatura¬ 
tion for an air bubble in water. The solid curve uses the satu¬ 
ration surface tension, 7 ^ = 0.072 N/m, and the dashed curve 
uses the supersaturation surface tension, 7 (s), with s* = 5. 


seen that the critical radius and the nucleation barrier 
are indeed reduced from their conventional values. The 
critical radius is nanoscopic for realistic parameter val¬ 
ues, and would move to even smaller values for larger 
reservoir supersaturation ratios, or if the spinodal super¬ 
saturation ratio had a smaller value. The critical radius 
is plotted as a function of supersaturation in Fig. [3] 

Explicit comparison can be made with measurements 
on nanobubbles. A supersaturation ratio of s = 3.6 and 
a spinodal supersaturation of s* = 5 gives a supersatu¬ 
ration surface tension of 7 (s) = 0.026 N/m and a critical 
radius of R CI ; t = 200 nrn. Alternatively, s = 4.4 gives 
7 (s) = 0.010 N/m and a critical radius of i? C rit = 60 nm. 
These are comparable to values obtained from measure¬ 
ments on nanobubbles, specifically from analysis of the 
difference between the nanoscopic and macroscopic con¬ 
tact angles^ with the former being obtained from from 
the tapping mode image profiles of nanobubbles In one 
series of experiments the average radius of curvature was 
R = 200 nm, and a surface tension of y(s) = 0.015 N/m 
was deduced from the contact angle. From another se¬ 
ries, R = 60 nm and y(s) = 0.019 N/m were obtained. 
An alternative interpretation of these results is discussed 
in the conclusion to the paper. 

From the curves in Fig. [2]it can be seen that the critical 
radius remains an entropy minimum when the supersat¬ 
urated surface tension is used. Hence the nanobubble 
remains unstable. In Appendix [C] the effect of using the 
local supersaturation (ie. the supersaturation in the solu¬ 
tion adjacent to the air-water interface, which is in diffu¬ 
sive equilibrium with the air inside the bubble), which is 
greater than the reservoir supersaturation, is discussed. 
Even in this case, and for a general non-linear model of 
the supersaturated surface tension, it is proven that the 
nanobubble is always unstable. Hence in general the su¬ 
persaturated surface tension alone cannot give a stable 
nanobubble. 

Despite the lack of stability, the decrease in surface 


tension with increasing supersaturation appears to solve 
one of the two major problems that nanobubbles pose for 
conventional thermodynamics. As has been previously 
argued^ it explains how the equilibrium radius comes 
to be of nanoscopic dimensions. 


IV. HEMISPHERICAL BUBBLE ON A 
SURFACE 


A. Mobile Hemispherical Bubble 


For a bubble on a solid surface, the optimum exterior 
contact angle is given by the Young equation 


7si = 7sg + Tig c°s( 7 T - 0 C ), or cos 6 C = ———. (4.1) 

Tig 

The liquid-vapor surface tension yi g was variously de¬ 
noted 7 , or 7 I, or 7 (s) above. For a hydrophobic surface, 
which is the present concern, the surface energies of the 
solid interfaces satisfy 7 sg < 7 si, or d c > 7r/2. Define 
A 7 = 7s g - 7 s i < 0. 

It should be noted that the difference in the solid sur¬ 
face energies is assumed fixed with a value determined 
from measurements on macroscopic water droplets on the 
given surface and the saturation surface tension. Con¬ 
sequently the optimum contact angle for nanobubbles, 
which is given by the Young equation, increases with in¬ 
creasing supersaturation ratio since the liquid-vapor sur¬ 
face tension decreases with increasing supersaturation ra¬ 
tio. This notion is the origin of the values of the surface 
tension given above, which were deduced from the mea¬ 
sured contact angle for nanobubbles 

Consider an air bubble with radius of curvature R, 
apex height above the surface of h, and contact rim radius 
r. The contact rim radius and apex height are related by 

r 2 +h 2 - 2 Rh = 0. (4.2) 


Hence r = \/2Rh — h 2 or h = R — \/ R 2 — r 2 . 

The interior contact angle is 9' = tv — 9, and the an¬ 
gle between the surface and the curvature radius at the 
contact rim is (f> = (n/2) — 6' = 9 — n/2. One has 
(R — h)/R = sin 4> = — cos 9, or cos 9 = —1 + h/R. 

The volume of the bubble is 


V = 


Rh 2 ~-h 3 


and the area of the liquid-vapor interface is 

Aw = 2t tR [l? - \/R 2 - r 2 1 = 2t xRh. 


(4.3) 


(4.4) 


The area of the solid-gas interface is of course A sg = 
7 r r 2 = 7 t (2 Rh — h 2 ). 

The total entropy is 


Stotai = S a (N,V,T) + - |U 

- 


(4.5) 
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Again one can use the ideal gas entropy for the bubble, 
S s (N,V,T) = 7Vfc B [ 1 ~ lnIVA 3 /U]. The total entropy 
is most conveniently written as a function of TV, R , and 
either h or r. 

Using Stota\(N, R, h\T), the derivatives are 

< 9 Stotal A* I. r \ 

~arT = — + r (4J3) 


dR 


and 


^^total 

dh 


El - L 
T T 


dV(R, h) 
dR 


7 i g dAi s (R, h ) Ay dA sg {R , h) 


T dR 

Ps _ P_ 

T T 


T 

2 7 ig + Ay 


7 xh — 


T 


dR 

27rh, 


(4.7) 


Ei- L 

T T 


dV(R, h ) 


yig 5Ai g (J?, h) Ay dA sg (R, h) 


T 


T dh 

E±_ P_ 

T T 

-^2nR-^2n[R-h\. 


dh 


7r [2 Rh - h 2 \ 


(4.8) 


If the total entropy is a maximum simultaneously with 
respect to curvature radius R and apex height h, then 
the left hand sides of Eqs (14.71) and (|4.8I) are zero. These 
two equations have simultaneous solution 

a R — h . ,. 

Ay = —y———=ycos 6, (4.9) 

R 

which is the contact angle condition, and 

Ps=P+ J[ ’ ( 4 -! 0 ) 
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FIG. 4: Total entropy of a mobile hemispherical air bubble in 
water on a hydrophobic surface, Ay = —0.036 N/m, following 
the mechanical equilibrium path. The solid curves are for a 
supersaturation ratio of s = 2 and the dashed curves are for 
s = 2.5. The lower curve of each pair uses the coexistence sur¬ 
face tension y+ = 0.072 N/m, and the upper curve of each pair 
uses the supersaturated surface tension y(2) = 0.054 N/m, 
and y(2.5) = 0.045 N/m. 


TABLE I: Geometry of a Mobile Hemispherical Bubble at 
the Critical Radius. The difference in solid energies is Ay = 
-0.036 N/m. 


s 

7 

-^crit 

T crit 

^crit 

9c 


N/m 

(/iUl) 

W 

W 

(deg.) 

2.0 

0.072 f 

1.44 

1.25 

0.720 

120 

2.0 

0.054* 

1.08 

0.805 

0.360 

132 

2.0 

0.048 # 

0.960 

0.635 

0.240 

139 

2.5 

0.072+ 

0.960 

0.831 

0.480 

120 

2.5 

0.045* 

0.600 

0.360 

0.120 

143 

2.5 

0.036 # 

0.480 

0 . 

0 . 

180 

3.0 

0.072+ 

0.720 

0.624 

0.360 

120 

© 

CO 

0.036* 

0.360 

0 . 

0 . 

180 

3.0 

O 

o 

to 

0.240 

0 . 

0 . 

> 180 


which is the Laplace-Young equation. These two results 
hold whether or not the entropy is maximal with respect 
to number (ie. whether or not Eq. m is zero). 

In Fig. HI the total entropy is plotted along the path of 
mechanical equilibrium, which is equivalent to imposing 
the Laplace-Young equation and the contact angle con¬ 
dition on the bubble as it grows. It can be seen that 
the total entropy is a minimum at the critical radius. 
The geometric parameters of a bubble at the critical ra¬ 
dius for these and other thermodynamic parameters are 
given in Table [I] An optimum contact angle of 180° or 
greater means that it is favorable for a planar gas layer 
to form between the solid surface and the liquid phase 
(drying transition). The results in Fig. [4] confirm that a 
hemispherical bubble mobile on a surface has the same 
stability characteristics as a free spherical bubble, which 
is to say that it remains unstable. 


4 Coexistence value. 

‘Supersaturated value using s+ = 5. 
* Supersaturated value using s* = 4. 


B. Pinned Hemispherical Bubble 

For the case of a hemispherical bubble on a planar sub¬ 
strate, if the contact rim of the bubble is mobile, (ie. as it 
is when it satisfies the contact angle condition) then its 
volume increases with increasing radius of curvature, just 
as for a free spherical bubble. In both cases the critical 
radius is unstable. However if the contact rim is pinned 
then an increase in the radius of curvature decreases the 
volume, and vice versa. In this case the critical radius 
becomes a point of stable equilibrium. 
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The physical origin of the stability is readily under¬ 
stood. Following closely the discussion of the instabil¬ 
ity of a free bubble given in the third last paragraph of 
the introduction, at the critical radius of curvature for a 
pinned bubble, a fluctuation to a larger radius of curva¬ 
ture along the path of mechanical equilibrium reduces the 
internal pressure and hence the gas chemical potential in¬ 
side the bubble. There is now a thermodynamic gradient 
that forces dissolved gas from the reservoir into the bub¬ 
ble, which increases the internal pressure. This creates 
a mechanical force imbalance, which acts to increase the 
volume of the bubble at the expense of the volume of the 
reservoir. So far this is the same as for a free bubble and 
for a mobile adhering bubble. But the difference for a 
pinned bubble is that an increase in volume reduces the 
radius of curvature, which counteracts the original fluc¬ 
tuation and restores the bubble to its original state. The 
converse occurs if the fluctuation is to a smaller radius of 
curvature. Hence the critical radius is a point of stable 
equilibrium for a pinned bubble. 

This physical argument that the pinned critical bubble 
is at an entropy maximum can be readily confirmed by 
numerical calculation. First one needs to derive the en¬ 
tropy derivatives, the vanishing of which give the critical 
radius. For a pinned (immobile) bubble attached to a 
hydrophobic surface, the rim radius r is constant. Hence 
needs to use r instead of h as an independent variable so 
that r can be held fixed during the derivatives with re¬ 
spect to R and N. Using the fact that h = R—\/R 2 — r 2 , 
the volume of the bubble may be written 
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FIG. 5: Total entropy of a pinned hemispherical bubble 
as a function of the radius of curvature for fixed contact 
rim radius r = 10 -7 m. The solid curves follow the me¬ 
chanical equilibrium path and dashed curves follow the dif¬ 
fusive equilibrium path. The lower pair of curves are for 
s = 2 and y(s) = 0.054 N/m, the middle pair use s = 3 
and y(s) = 0.036 N/m, and the upper pair use s = 4 and 
7 (s) = 0.018N/m (s* = 5 in all cases). The arrows denote 
the entropy maxima, which are at the critical radii. 


and 

dS total (N,R\r,T) 

dR 

-p s p-i dV(R,r) 7 i g dA\ g (R, r) 

.T T\ dR T dR ' 1 ' ’ 


V = 


7r 


7r 


Rh 2 - -h 3 
3 

-R 3 - -R 2 yjR 2 -r 2 
3 3 


(4.11) 



and the area of the liquid-vapor interface is 


A lg 


2ttR 


r-^/r 2 



(4.12) 


Hence at the bubble stationary point when the radius 
derivative vanishes (mechanical equilibrium for a pinned 
contact rim) one has 


Ps 


dAi g (R, r) 3R 
P 7g dR dV{R,r) 

4 VR 2 - r 2 -4 R + 2 r 2 /R 
Pr + 7lg 2R\JR 2 - r 2 - 2 R 2 + r 2 ' 


(4.18) 


The area of the solid-gas interface is of course H sg = irr 2 . 
The derivatives are 

= *t 2 f+ ’ J l fi 

dR 3 3\/ R 2 - r 2 

(4.13) 


dA lg (R,r) 

dR 


= 47 tR — 2i r\/ R 2 — r 2 — 


2t tR 2 


y/R 2 ^' 


f, (4-14) 


dA sg (R,r) 

dR 


= 0. 


(4.15) 


With these, the derivatives of the total entropy at fixed 
r are 


dS total (N,R\r,T) _ -p s p 
dN T + T’ 


(4.16) 


If R ss r, then p s ~ p T + 2 71 g /R, which is the usual 
Laplace-Young expression. More generally, the left hand 
side is p s = Nk^T/V(R,r). Inserting this into the above 
and taking the volume over to the right hand side gives 
an explicit equation for N(R) along the path of mechani¬ 
cal equilibrium. This may be inserted into the expression 
for the total entropy and the latter plotted as a function 
of the radius of curvature for R> r. For the case of diffu¬ 
sive equilibrium, N = sp^V(R,r). Either the saturated 
surface tension 7 I or the supersaturated surface tension 
7 ( 5 ) may be used. 

Figure O shows the total entropy as a function of the 
radius of curvature when the contact rim is pinned at 
r = 1CT 7 m. Unlike the previous cases, the total entropy 
is now a maximum at the critical radius. The maximum 
is rather broad and extends down almost to the pinned 
radius itself. For supersaturations of s = 3 and s = 4, 
the total entropy maximum is positive, which indicates 
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FIG. 6: Total entropy of a pinned hemispherical bubble (s = 
2, y(s) = 0.054 N/m, r = 10 _7 m) in the neighborhood of 
the critical radius of curvature and density (±20% in each 
dimension). 


that it is favorable to have a pinned bubble compared to 
no bubble at all. 

Figure O is a surface plot of the total entropy for the 
pinned hemispherical bubble over the curvature radius- 
density plane. It can clearly be seen that the critical 
radius and density is a point of global entropy maximum. 
In other words the pinned (immobile) bubble is stable. 


V. CONCLUSION 


In the introduction four contradictions between classi¬ 
cal nucleation theory and nanobubbles were pointed out: 
(1) in the classical theory bubbles are unstable whereas 
nanobubbles appear stable, (2) the classical critical ra¬ 
dius is 1-2 orders of magnitude larger than the nanobub¬ 
ble radius of curvature, (3) the classical internal pressure 
of nanobubbles is 1-2 orders of magnitude larger than 
atmospheric pressure, and (4) the macroscopic contact 
angle is significantly larger than the measured nanobub¬ 
ble contact angle. 

It was shown here that two non-classical ideas are 
necessary and sufficient to account for the behavior of 
nanobubbles, namely the dependence of surface tension 
on supersaturation, and the pinning of the contact line 
of the nanobubble on the surface. 

It has previously been pointed out that bubbles can 
only be in diffusive equilibrium with a supersaturated 
solution^ and that the supersaturated surface tension is 
smaller than the saturated surface tension^ These facts 
reduce the size of the critical radius and the magnitude 
of the internal pressure from those predicted by classical 
nucleation theory. This largely resolves the second and 
third problems enumerated above since, depending upon 
the parameters, it is possible to get quantitative agree¬ 
ment with the measured radius of curvature of nanobub¬ 
bles and the deduced value of the surface tension^ (see 
below) for apparently reasonable values of the parame¬ 


ters (actual supersaturation ratio of the solution and the 
spinodal supersaturation ratio). The origin of the super¬ 
saturation has been variously attributed to the heating 
of the solution and the surface by the laser diode, the 
piezo-electric crystal, and exothermic mixing during the 
exchange of ethanol for water, and also to the entrain¬ 
ment of air by the surfaces^ In any case, whatever the 
origin of the air, once stable nanobubbles are present, the 
inescapable thermodynamic conclusion is that the solu¬ 
tion is supersaturated. 

The new contributions here on this point were the ex¬ 
plicit calculations for free and adsorbed bubbles, and the 
mathematical proof that the state-dependent surface ten¬ 
sion on its own can never give a stable nanobubble, Ap¬ 
pendix 1C 21 

On the first point concerning stability, it was confirmed 
that neither a free spherical bubble nor a mobile hemi¬ 
spherical bubble on a solid surface were stable. However, 
if the bubble was pinned with fixed contact rim, then 
it became stable. The physical reason for the change 
in stability is that whereas both a free spherical bubble 
and a mobile hemispherical bubble have a volume that 
increases with increasing radius, the volume of a pinned 
hemispherical bubble decreases with increasing radius of 
curvature. This is sufficient to confer stability. 

The fourth point is resolved by either or both of these 
non-classical notions. The decrease in surface tension 
from the saturated value of a macroscopic droplet will 
increase the contact angle provided that the solid surface 
energies remain unchanged. This observation was the ori¬ 
gin of the estimate of the surface tension of nanobubbles 
given in Ref. [l8i . and it appears to be consistent with their 
measured radii of curvature. However, if the nanobubble 
is pinned, then the contact angle will also increase as the 
radius of curvature grows toward its critical value. 

For example, for the case of a macroscopic contact an¬ 
gle of 120°, (Ay = —0.036N/m), a mobile bubble using 
the saturated surface tension yt = 0.072 N/m has the 
same contact angle at the critical radius. At s = 2 
and y(s) = 0.054 N/m, a mobile bubble has contact 
angle 132° at the critical radius. A bubble pinned at 
r = 100 nm has contact angle 176° for yt = 0.072 N/m, 
and 172° for y(s) = 0.054 N/m, each at the respective 
critical radius (entropy maximum). 

Pinning evidently has a more dramatic effect on the 
contact angle than supersaturation. The fact that the 
measured contact angles for nanobubbles are larger than 
that of macroscopic water droplets on the same surface 
supports both non-classical notions and on its own can¬ 
not decide between them. Conversely, the calculations in 
Ref. [I|| of the surface tension of nanobubbles from the 
measured contact angles are probably not quantitatively 
reliable because they neglect the possibility of pinning. 

Two questions now arise: is pinning the only mecha¬ 
nism that will confer stability to a bubble? And what is 
the physical origin of pinning? 

It has previously been shown that a finite one- 
component system leads to stable bubbles^ and it is 











FIG. 7: Sketch of a nanobubble pinned by its deformation of 
the solid surface (not to scale). 

of interest whether similar results hold for nanobubbles 
derived from air dissolved in water. A model of a finite 
reservoir based upon the limited diffusion of air in water 
in investigated in Appendix [D] This model does yield an 
entropy maximum and stable bubbles at radii beyond the 
classical critical radius, but only on unrealistically short 
time scales. It is concluded that on experimental time 
scales the reservoir of dissolved gas is effectively infinite 
due to diffusion. Hence it appears that pinning is the 
only viable source of stability for nanobubbles. 

The most obvious sources of pinning are permanent 
heterogeneities on the surface. However, most measure¬ 
ments on nanobubbles are performed on surfaces deliber¬ 
ately chosen to be smooth and chemically homogeneous. 
This does not rule out random heterogeneities in any 
one case, particularly when isolated nanobubbles are ob¬ 
served. Because nanobubbles are so small, a very low 
surface density of heterogeneities may be sufficient to pin 
the contact line at a few points, which would immobilize 
the whole rim. Conversely, pinning at a few nanoscopic 
heterogeneities likely has negligible effect on macroscopic 
bubbles or droplets. 

Another possible source of pinning is the electric dou¬ 
ble layer repulsion between neighboring nanobubbles. 
The images of nanobubbles obtained by Tyrrell and 
AttardM show them to be close packed on the surface, 
exactly as one would expect for the surface nucleation 
from a supersaturated solution. Further, the contact 
lines are often irregular and non-circular, whereas a mo¬ 


bile bubble would be expected to have a circular con¬ 
tact line. The nanobubbles are negatively charged and 
their morphology changes with pH. Although the electric 
double layer repulsion between neighboring nanobubbles 
would prevent their contact lines from growing, it would 
not prevent them from shrinking. One concludes that it 
is probably not a general a source of stability. 

Even in the absence of permanent heterogeneities or 
neighbor interactions, it is possible that self-induced pin¬ 
ning of the contact line may occur due to the elastic or 
viscoelastic deformation of the substrate by the nanobub¬ 
ble, as is sketched in Fig. [7J The deformation is due to 
the twin effects of the excess pressure inside the nanobub¬ 
ble and the tension of the air-water interface applied at 
the contact line. Depending upon the relaxation time of 
the substrate, the barriers to relaxation, and the driving 
forces for horizontal motion, the contact rim may be effec¬ 
tively pinned at the deformation. It is clear that the cost 
of deformation of the base and rim of the substrate favor 
a small contact radius, whereas the cost of curvature of 
the contact rim favors a large contact radius. The latter 
effect is much more significant for a nanobubble than for 
a macroscopic bubble. 

In conclusion, this paper has shown that two non- 
classical notions are necessary and sufficient to account 
for the size and stability of nanobubbles: the surface ten¬ 
sion depends upon the level of supersaturation of the 
solution, and pinning of the contact rim confers thermo¬ 
dynamic stability. 

These appear to solve the theoretical thermodynamic 
problems posed by nanobubbles. They indicate the direc¬ 
tion for more experimental work to further characterize 
nanobubbles and to establish that control over the phe¬ 
nomena that is a prerequisite for industrial or technolog¬ 
ical exploitation. Quantitative measurements with con¬ 
trolled supersaturation, either by means of de-aeration, 
over-gassing, or temperature control, represent one line 
of enquiry. Measurements with substrates of different 
deformability, roughness, or polymer extensivity are also 
suggested by the present results. 
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Appendix A: Two-Component Formulation 


At the more fundamental level, a bubble and sur¬ 
rounding liquid should be treated as two coexisting two- 
component mixtures, one rich in air, and the other rich 
in water. Let air be denoted by the subscript 1 and water 
by the subscript 2. The total entropy is 


Stot{N 1 ,N 2 , V\ni,p 2 ,p,T) 

= S s (N u N 2 ,T) + ^ + ^ 


pV 

~T 


(A.D 


Equilibrium is characterized by the vanishing of the 
derivatives, namely 


Msl ~ Ml) Ms2 — M2 

(A.2) 

and 


27 

p - = t+ r- 

(A.3) 

At saturation one has 


p = p\ +p\, 

(A.4) 


where p is atmospheric pressure. For a supersaturation 
ratio of s, the liquid solution phase (ie. the reservoir) is 
in equilibrium with a gas phase of partial air pressure sp\ 
and partial water vapor pressure of p \. When the bub¬ 
ble is in diffusive equilibrium with the solution, the sum 
of these must give the internal pressure of the bubble. 
Hence the critical radius is given by 


27 

sp\ +p\~p 
27 

(s - 1 )p\ 


(A.5) 


The difference between this and the one-component for¬ 
mulation is that the partial vapor pressure of air at sat¬ 
uration, p\ ft 0.8 p, has been replaced by atmospheric 
pressure, p. 


Appendix B: Unstable Critical Radius 

As was shown explicitly in Fig. |T] above, the critical 
radius is a saddle point, which is an unstable equilibrium 
point. This may be proven mathematically by showing 
that a bubble does not satisfy the general stability re¬ 
quirements. 


1. General Stability Condition 

Since the probability of a fluctuation in number or vol¬ 
ume is proportional to the exponential of the entropy, 
two conditions necessary for stability are that the corre¬ 
sponding second derivatives of the entropy are negative, 

Snn < 0, and Syy < 0, (B.l) 

at least at a stable stationary point. More generally, 
the necessary and sufficient condition for stability is that 
both eigenvalues of the Hessian matrix of second deriva¬ 
tives are negative. The eigenvalues are determined from 
the characteristic equation, 

Snn — A Snv 
S yN Syy — A 

= [Snn — A] [Syy — A] — Sy N 
= A 2 — A [Snn + Syy] + [SawS yy — S^y] (B.2) 

whose roots are 

A± = - | [Snn + Syy] (B.3) 



The argument of the square root is always non-negative. 
The only way for both eigenvalues to be negative is if 
the square root is smaller than the magnitude of the first 
term, which occurs when the second bracketed term in 
the square root is positive. Hence in addition to the 
negativity of the two pure second derivatives, stability 
is ensured by the positivity of the determinant of the 
Hessian matrix, 

SnnSw — S 2 nv > 0. (B.4) 


2. Bubble Instability 

Specifically for a bubble, the second derivatives of the 
total entropy are 


Snn 


—fee 


Syy 


-pk B 

V 

~Ps , 


2j 1 

TR? AnR 2 
27 


(B.5) 


(B.6) 


VT 3 VTR 
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and 


The derivative with respect to radius is 


Snv = -y. (B.7) 

Clearly Snn is always negative. At the equilibrium 
point, the Laplace-Young equation shows that Syv = 
— (Pr/VT ) — (2y /VTR) + (2y/3 VTR), which is negative. 
These two conditions are necessary but not sufficient con¬ 
ditions for the equilibrium point to be stable. 

The determinant of the Hessian matrix is 


cib'total 

d R 


'Ps _ Pr] dV{R) 

.T T\ dR 

7 dA(R) i{p)Adp(R,N) 
~ T dR T dR 


'Ps 

.T 


_ Ih 
T 
7 (P) 
T 


AnR 


8ttR + 7 r (p) 


12nRp 

T 


(C.3) 


SnnSvv — -Sivy 


M. / Zi L p + _J?2_ 

V 2 \ P [ H 3 k B TR 


-fcj 27 

V 2 3pk B TR' 



which is always negative. As shown above, it is necessary 
for stability that this determinant be positive. 

One concludes from this that the bubble’s stationary 
point, the critical radius R c rit and density p C rit, is a saddle 
point, which is a point of unstable equilibrium. 


Appendix C: Local Supersaturation 


This vanishes when 


P s = Pr + 


2 b(p) 

R 


3pY(p) 

R 


(C.4) 


Like the usual Laplace-Young equation, this says that the 
internal pressure of the air bubble is larger than the exter¬ 
nal pressure of the liquid, since q(p) > 0 , and 7 '(p) < 0 . 
Since the usual Laplace-Young equation uses the coexis¬ 
tence value for the surface tension, 7 ^ = 7 (/c 4 ) < 7 (/?), 
whether the internal pressure predicted by the present 
theory is larger or smaller than that predicted by the 
constant surface tension theory depends on the specific 
surface tension function. 


1. Critical Radius for Local Supersaturation 


The simulations that justify the formula for the super¬ 
saturated surface tension were carried out with diffusive 
equilibrium established across the interface. In homoge¬ 
neous nucleation theory it is more traditional to plot the 
total entropy (equivalently free energy) as a function of 
radius along the path of mechanical equilibrium. In this 
case diffusive equilibrium does not hold except at the crit¬ 
ical radius. One could argue that there is always a local 
diffusive equilibrium in the immediate locality of the bub¬ 
ble air-water interface, and that this is what determines 
the value of the surface tension. In this case one would 
use 7(siocai) = 7(p)> where si oca i(p) = p/p\ = Ps/p|, for 
the bubble on the path of mechanical equilibrium. 

If one uses this bubble-dependent model for the super¬ 
saturated surface tension, the total entropy is 


S total {N,R,T ) 


S s (N,V,T) + ^ 

_ l(p)A 
T ’ 


PiV 

T 

(C.l) 


2. Instability of the Critical Radius for Local 
Supersaturation 

At the bubble critical radius (p C rit> ^?crit), both deriva¬ 
tives vanish, Sjv(p C rit, -Rcrit) = 0 and -SV(/9 C rit, -Rcrit) = 0 . 
Individually they give a non-linear function for the opti¬ 
mum density 


p = sp f e- 3 ' 1 ' M/kBTR , 


and 


t , 2 7 (p) - 3 fn'{p) 
P = P + - k^TR - 


(C.5) 


(C. 6 ) 


The simultaneous solution of these gives the critical ra¬ 
dius and density. If the derivative of the surface tension 
is negative, 7 '(p) < 0 , then the first of these says that the 
gas density inside the bubble is greater than the super¬ 
saturation vapor density. The second says that the gas 
density inside the bubble is greater than the coexistence 
density, both terms of the correction being positive. 

Solving instead for the critical radius one obtains 


with p = N/V and p r = p^bT. For full generality, and 
to make an important point, it will not here be assumed 
that 7 (p) is a linear function of density. 

The number derivative is 


•AS'total 

dN 


-Ps br 

T + T 


7 '{p) 


A 

VT' 


(C.2) 


R = 


—37 '{p) 
k B T 



— = 27 (p) - 3pj{p) 
k B T[p - pt] 


(C.7) 


(C. 8 ) 


When this derivative vanishes, the chemical potential in- The second of these gives R(p) for p £ along the 

side the bubble is higher than that outside. path of mechanical equilibrium. 
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Since the critical radius must vanish at the spinodal, 
R(Pg) = 0 > this says that 

tO4) = \p\i{p\)- ( c - 9 ) 

It is not possible to satisfy this in the linear model, since 
the left hand side is positive or zero and the right hand 
side is strictly negative. Hence the local supersaturation- 
dependent surface tension must be a non-linear function 
of the internal gas density. Since in general Y(Pg) < 0, 
and this result says that Y(p\) > 0 , one must have that 

7 "(P) > 0, (C.10) 

assuming that the interval is relatively small. 

This result holds only for the local supersaturation- 
dependent surface tension. It will now be shown that this 
result precludes the critical radius from being a point of 
stable equilibrium. 

The second derivatives of the total entropy are 


_ -ks _ 3y "(p) 

N VTR ’ 


(C.ll) 


The determinant should be evaluated at the bubble 
critical point, which is denoted by an over-line. For typo¬ 
graphical simplicity, the density argument of the surface 
tension will be suppressed. One has 


V 2 

-Tt-SnnSw = 

kl 


and 


V 2 

k B 


&NV ~ U + 


-i 3 7 " ] 
p Rk B T J 

2j — 6pY — 9p 2 Y' 


x <; -p+ 
1 - 


3Rk B T 
27 — 6pY — 9p 2 Y' 3 p 7 ‘ 


3RpkYT Rk B T 

277 " - 6pYY' - 9p 2 Y' 2 


7 


(. Rk B Tf 


3/07" y 


(C.18) 


(C.19) 


= 1 + 


Rk B T Rk B T J 
2Y + QpY' 6 pYY' + Y 2 + $p 2 Y' 2 


Rk B T 


(. Rk B Tf 


For a stable bubble critical point, the determinant should 
be positive 


~P k b 27 (p) - 3pY{p) 1 

V TR 2 4t tR 2 

_ Y(p) + 3pY'{p) p 

TR V 

-pk B , 2 j(p) - 6pY(p) - Vp 2 Y'(p)rr, 

“f” + WTR 

and 


^ — 2 "/ + 6pY 277" — 6pYY ' ~~ 9p 2 Y' 2 

3Rpk B T (Rk B T) 2 

27 ' + 6 / 97 " 6p7 , 7 ,/ + 1' 2 + 9 p 2 7" 2 

Rk B T ( Rk B T ) 2 

-27 - 18p 2 7" 277" + Y 2 

3 Rpk B T (Rk B T) 2 ' 

This gives 


Svv = 


k B 37 '(p) 1 37 " (p) p 

Snv = v + -TWYYtf + -tTv- (ai3) 

Clearly Snn is negative if 


„ —k B TR 
7 (P) > — 

3 P 


(C.14) 


Although this bound could be satisfied by Y'{p) = 0 , this 
possibility is ruled out by another bound below. 

Also Svv is negative if 

0 > -n B T + (C . 16 ) 


3 R 


which can be written as 


. ^ -3 Rpk B T + 2 7 (p) - 6pY(p) , c ^ 

7 ip) >-—}- 1 (C.16) 


or as 


7 < ( Rk B T) 2 


- 27 - 18 p 2 7 " 


277" 


3 Rpk B T ( Rk B T) 2 \ 


(C.21) 


This has a solution if, and only if, the bracketed term is 
positive, which means that 


-2 7 


> J 7 '' 


— 2 —// 


18/9 7 ' 


3Rpk B T ( Rk B T ) 2 3Rpk B T 

27 + 6pRk B T_„ 

~ (Rk B T) 2 1 ’ 


or 


7 < 


-27 ( Rk B T ) 2 

3pRk B T 27 + PPpRkrQT 
— Rk B T /3 ~p 
1 + 3pRk B T/j 


(C.22) 


(C.23) 


Combining this with the condition found above for 
Snn < 0, 7 " > — Rk B T/3p } gives 


Y(p) > ~ 3R P kBT + ~ ®p 2 Y'(p) 17 ^ 

6p 


RkB R 3p < -i> < (C.24) 

1 + 3pRk B T/*? 3p 
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This rules out the possibility that 7 " = 0, which is to 
say that the local supersaturated surface tension cannot 
be a linear function of density if a bubble is to be stable. 
This agrees with the conclusion already made above on 
different grounds. 

This result shows that the second derivative of the sur¬ 
face tension has to be negative if the free bubble is to 
be stable. But Eq. (1C.101) showed that the second den¬ 
sity derivative of the surface tension had to be positive 
7 "(p) > 0. These two results prove that it is not possi¬ 
ble to confer thermodynamic stability on a bubble at the 
critical radius by invoking a local supersaturation surface 
tension. 



R ( m ) 


Appendix D: Finite Reservoir 


In order to model the finite diffusion of air in water, one 
can consider two reservoirs: an infinite volume reservoir 
of atmospheric pressure p, and a reservoir of fixed volume 
V r . The total volume Vtot = V + V r is variable and can 
exchange with the infinite volume pressure reservoir. The 
total number of of air molecules is fixed, N tot = N+N r = 
spV r /K , where I\ = 2.3 x 10 -19 J is Henry’s constant for 
air in water. The idea is that due to limited diffusion over 
a specified time scale, only air dissolved in the water of 
the fixed volume reservoir can exchange with the bubble. 

One can define the fixed reservoir volume in terms of 
a diffusion length 

k = 1d ~^ (D- 1 ) 

where D = 2x 10 -9 m 2 /s is the diffusion constant for air 
and water, and t is some experimental time scale. 

The spherical bubble or sub-system is characterized 
by number N, volume V = 4ttR 3 /3, and area A = 47 tI? 2 . 
With these the total entropy is 

Stotal = Ss(N, V, T) + S t (N t , V t ,T) - j;A 

- |;Vtot + const. (D.2) 

The constant is chosen to make the entropy vanish at 
zero radius. 

As usual the sub-system (bubble) is an ideal gas 


S s (N,V,T) = k B N 


1 — In 


N A 3 ' 


(D.3) 


The fixed volume reservoir containing the dissolved air 
can be treated as an ideal gas in a uniform external field 


S T (N r , Vr, T) = k B N T 


1 — In 


N r A 3 ' 

”TT. 


N t £ 


(D.4) 


Here e is the molecular solvation energy. The effective 
chemical potential that corresponds to this is 

/A = £ + /c B Tln p S oiA 3 , (D.5) 


FIG. 8: Total entropy of a spherical bubble in a finite reser¬ 
voir on the path of mechanical equilibrium (s = 2, 7* = 
0.072 N/m). The solid curve is for an infinite reservoir, the 
dashed curve is for a diffusion time t = 0.035 s, and the dotted 
curve is for t = 0.01 s. 


where p so i = N T /V T is the concentration of dissolved gas. 
Henry’s law may be written for the concentration of dis¬ 
solved gas in the form 

p sol = K~ x p = K- l k B TA~ s e^ k * T , (D. 6 ) 


or 


li = ksT In 


K P sol 

k B T 


(D.7) 


Recall that the saturation air vapor pressure has been 
taken to equal the external atmospheric pressure. Equat¬ 
ing these two expressions for the chemical potential yields 

£ = k B T hi^-- (D. 8 ) 


Because the reservoir volume V r is fixed and indepen¬ 
dent of the bubble, the derivative of the total entropy 
with respect to radius vanishes when the Laplace-Young 
equation is satisfied, 

P s =Pext + (D.9) 

This gives the path of mechanical equilibrium. 

Figure [S] plots the entropy using this model for a fi¬ 
nite, diffusion-linked reservoir. The results are typical of 
the model. For small bubbles the reservoir is effectively 
infinite and all the curves coincide. For larger bubbles 
and small reservoirs (ie. short experimental time scales), 
the total entropy monotonically decreases and there is 
no critical radius. For larger reservoirs (ie. longer experi¬ 
mental time scales), there is a critical radius and beyond 
this a stable radius where the entropy is a maximum. 
The decrease in entropy beyond this maximum is due to 
the depletion of the reservoir. Obviously, the radius of 
maximum entropy will increase with increasing reservoir 
size. It should be noted that the typical lifetime observed 
for nanobubbles is on the order of 10 3 - 10 4 s. 














